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Abstract. We investigate the usage of highly efficient error correcting codes of 
multilevel systems to protect encoded quantum information from erasure errors and 
implementation to repetitively correct these errors. Our scheme makes use of quantum 
polynomial codes to encode quantum information and generalizes teleportation based 
error correction for multilevel systems to correct photon losses and operation errors 
in a fault-tolerant manner. We discuss the application of quantum polynomial codes 
to one-way quantum repeaters. For various types of operation errors, we identify 
different parameter regions where quantum polynomial codes can achieve a superior 
performance compared to qubit based quantum parity codes. 


1. Introduction 


A quantum erasure channel replaces a qubit (qudit) with an “erasure state” that is 
orthogonal to all the basis states of a qubit (qudit) with a certain probability, thereby 
erasing the qubit (qudit) and enabling the receiver know that it has been erased [I]. 
Physically, erasure errors may occur in various situations, such as leakage to other states 
[2| [3 [4], atom losses [5], and photon losses [6} [7 [8] [9]. For ion-trap systems, leakage 
processes occur when the qubit moves out of the idealized two-level sub-space to a 
larger space [4]. For quantum memories with optical lattices, back ground gas 
collisions can eject the atoms leading to atom losses Photon losses occur in linear 
optical quantum computing schemes [6] due to absorption in optical interconnects or 
optical fiber. Undoubtedly, protecting quantum information from erasure errors is a 
significant challenge for practical quantum computation and long distance quantum 
communication. 

Specifically, for long distance quantum communiction through optical fibers, 
photon losses lead to an exponential penality in resources and time. The exponential 
penality can be overcome by establishing intermediate repeater stations and actively 
correcting for erasure and operation errors at these stations. Three generations of 
quantum repeaters have been proposed based on the different approaches used to 
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correct erasure and operation errors [10]. The first generation employs heralded 
entanglement generation between neighboring repeater stations to correct erasure errors 
and entanglement purification to correct operation errors {11}. Heralded entanglement 
generation needs two-way classical communication between neighboring repeater 
stations, while entanglement purification needs two-way classical communication 
between remote repeater stations. The second generation employs 
heralded entanglement generation to correct erasure errors and quantum error correction 
to correct operation errors. Quantum error correction doesn’t require any form of 
two-way classical communication. The third generation uses quantum error correction 
to correct both loss and operation errors, and avoids any form of two-way classical 
communication between repeater stations, thereby rendering ultrafast communication 
over transcontinental distances [16} [1’7) [18} [79] [20] (21). Since erasure errors are actively 
corrected in these repeater schemes, it is crucial to investigate quantum error correcting 
codes that can correct erasure errors very efficiently 23} (24) [25]. So, far only quantum 
parity codes (QPC) have been optimized for third generation quantum repeaters [17]. 

Due to the quantum no-cloning theorem [I], no error correcting code can correct 
erasure errors deterministically when the erasure rate is above 50%. There have been 
significant advances in searching for quantum codes that can correct up to 50% erasure 
error rates. Varnava et. al. [22] showed that by using tree-like cluster states for 
encoding, one can correct erasure errors when erasure rate is close to the 50% bound for 
one-way quantum computation. Stace and Barrett [24| demonstrated that surface 
codes can also correct erasure errors when error rates are close to 50%. However, qubit 
based codes often require large code size to enable the correction of a large fraction of 
erasure errors. 

Quantum error correcting codes of higher-dimensional systems provide a promising 
alternative to qubit based encoding schemes to correct erasure errors. For example, 
quantum polynomial codes (QPyC) are a class of CSS codes that were introduced 
in the context of fault-tolerant quantum computation and shown to be useful for 
constructing threshold quantum secret sharing schemes [28]. One can encode a secret 
qudit into 2k +1 qudits (with prime dimension d > 2k + 1) and distribute one qudit 
to each of the many parties, so that at least (k +1) of them should get together to 
reconstruct the secret. This makes the [[2k + 1,1, + 1]|a QPyC a good choice for the 
correction of erasure errors up to a fraction of k/(2k + 1) + 50% for a large k. 

This paper is organized as follows: We first perform a comparison between |[3, 1, 2]]3 
code and the [[4, 2, 2]] code. In section II, we investigate the ability of general QPyC 
to correct erasure errors and compare them to surface codes. In section IV, we show 
that QPyC can be used for third generation (or one-way) quantum repeaters. We 
also compare our quantum repeater scheme with other schemes based on quantum 
parity codes in the presence of operation errors. Here, we identify the 
parameter regimes where QPyC performs better than QPC. The price for this improved 
performance against erasure errors is that more complex multi-mode operations must 
be implemented for encoding and readout operations. In section V, we discuss the 
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key experimental techniques needed for the physical implementation of our scheme and 
provide potential experimental procedures based on atom mediated photonic gates for 
multilevel systems. 


2. Three-qutrit code vs Four-qubit code 


To illustrate how the error correction for erasure errors works, consider the four-photon 
[[4, 2, 2]] code that maps two qubits into the logical states [25] [29], 


00), = | 00) + ]11))12({00) + |11))34, 
01) = 5(\00) — [11))12(100) —|11))se (1) 
10)1, = 5(J01) + [10))12(101) + |10))s4 
11), = 5(J01) ~ [10))12(101) —[10))s 


where the subscripts 1-4 denote qubits 1-4. Here, |0) and |1) corresponds to a single 
photon occupying two different modes such as polarization or time-bin states ’01’ and 
10’ respectively. The [[4,2,2]] code requires a total of four pairs of modes (eight modes) 
and four photons. The loss of a photon leads to the vacuum state of the associated 
pair of modes, which corresponds to an erasure error. Suppose we transmit the encoded 
state through a channel and it undergoes one photon loss (i. e. erasure of one of qubits), 
then the quantum error correcting code enables the reconstruction of the encoded state 
as follows. First, we use quantum non-demolition measurement to extract the total 
excitation number for each pair of modes (without destroying the qubit), which will be 
one if there is no photon loss, or zero if the photon is lost. Hence, the photon loss can 
be identified as an erasure error of the associated qubit. For example, if the first qubit 
is erased, to reconstruct the encoded state we apply two CNOT gates betweeen qubits 
3 and 2, and qubits 4 and 2 respectively. Then, we measure qubit 2 in the Z basis and 
and use the measurement outcome to reconstruct the logical state. This can be seen 
by studying the logical operators of the [[4, 2, 2]] code, x =I XIX, x? = 1ZTZ. 
ZY) = IITZZ and Ze) = IIXX. After the CNOT gates, the logical operators are 
transformed into IJIX, IZZI, IIZZ and IIXX respectively. After the gates, a Z- 
measurement on the second qubit is needed to decode the [[4, 2, 2]] code. We refer the 
readers to Ref. [25], where an alternative non-destructive method has been proposed 
to recover the logical qubits. It has been proven that we need at least four qubits to 
correct one erasure error [30]. 

Alternatively, one can also correct an erasure error using a [[3, 1, 2]]3; code which 
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encodes a logical qutrit into three physical qutrits as [28], 


1 

0), = Wai 000) + |111) + |222))193, (2) 
1 

Dee We 012) + |120) + |201))193, 
1 


where the subscripts 1-3 represent qutrits 1-3. Here, we note that each qutrit represents 
different time bin states and not photon number. More specifically, |0), |1) and |2) 
represent temporal modes ’001’, ’010’ and ’100’ respectively. So, the [[3,1,2]]3 code 
requires a total of three triplets of modes (nine modes) and three photons. Photon loss 
leads to the vacuum state of the corresponding triplet of modes, which can identified as 
the erasure error of the associated qutrit. If qutrit is erased, it is possible to reconstruct 
the encoded qutrit by performing an addition modulo 3 operation between the other 
two qutrits. For example, for the incoming state a|0);+|1)r+/|2)z, if the first qutrit 
is erased, then the state of second and third qutrits is given by 


P23 = S- \Widos (Wal, (3) 


1=0 


where, 


[eo) = a|00)23 + B]12)03 + 721) 23 
|v) = a|11)23 + B]20)23 + 702) 23 
|t2) = a|22)23 + B|01)23 + 7|10)23. 


Using two addition modulo 3 operations{] between qutrits 2 and 3, and then between 
qutrits 3 and 2, we can reconstruct the encoded state. 

We now compare the |/3,1,2]]; code and the [[4, 2, 2]] code for different loss rate 
per photon, p;. Since each logical qutrit can carry log,3 qubits of information, we 
multiply a factor of log.3 when computing bits/photon and bits/mode for [[3, 1, 2]]3 
code. As illustrated in Fig. the bits/photon is always higher for the [[3, 1, 2]]s 
code compared to the [[4, 2, 2]] code for all loss rates py. With regard to the bits/mode, 
for large loss rates (p, > 42%), the [[3,1,2]]; code performs better than the [[4, 2, 2]| 
code. In principle, by concatenating the |[3, 1,2]]3 code, it is possible to supress more 
erasure errors. However, to be resource efficient we will consider the generalization of 
the [[3, 1, 2]]; codes - QPyC in the forthcoming section. 


3. Quantum Polynomial codes 


The qutrit [[3,1,2]]; code can be generalized to d-level (d is a prime number) qudit 
system as a |[2k+1,1,k+1]]q code, which encodes one logical qudit into 2k +1 physical 


t An addition modulo 3 operation refers to a SUM gate here is described in detail in the section IV. 
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Figure 1. a) A comparison between bits/photon that can be achieved with [[3, 1, 2]]s 
code and [[4, 2, 2]] code respectively. b) A comparison between bits/mode that can be 
achieved with [[3, 1,2]]3 code and [[4, 2, 2]] code respectively. 


qudits and can correct up to k erasure errors [28]. The QPyC code is a CSS code which 
can be obtained from classical [2k + 1,1,4-+ 1] polynomial code. The encoded states of 
the QPyC code are given by [28], 

[s)z = [p(C; 25); D(C, 21) a4 D(C; L2%))- (4) 


cp=s,|ceF*t! 


Here, s labels the logical states with s € {0,1,2...d—1}. The labels xo,24,...,72, are 2k+1 
distinct elements of F, where F = Z. d is the first prime number that is no less than 
2k +1. The polynomial is given by p(c,t) = 2 cyt? with ¢ = (co, c1,...ce) € F**?. 
For k = 1, d= 8, it can be seen that 


2 
ls), = Fs lco)|co + €1)|co + 2c1), 


co=0 


(5) 


yielding the [[3, 1,2]; code (see Eq. (3p). The success probability of recovering the 
encoded quantum information after it has undergone erasure errors is given by 


k 


2k+1 _3 


j=0 


pebyc 


SUCCESS 


(6) 


Note that for a reasonable code size, it is possible to achieve substantially low (e. g. 
10~°) failure rate Pray = 1 — p@PyC even in the presence of high erasure rate pj. Note 


that with about 40 qudits one can suppress the failure rate to 10~° for 20% erasure rate. 
It can be seen from Eq. (6) that at p, = 50%, P2EYS = 1/2 independent of the code 


success 
size. Further, note that the success probability of error correction has a phase transition 


behaviour with k + 00, P2PYS — 1 for p, < 50% and P2PYS — 0 for p, > 50%. One 


success success 
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can obtain the critical exponent of phase transition by noting that for k — oo, around 
the fixed point PQPYS (k, 4) = x, the asymptotic expression is 


success \'¥) 9 


POPYC (ki, py) = 1 - ( oe ) -an—1 VAT (k + 2) 


ar(k + 3/2) 
op, 4k2Fi (2,1 — kk +3;—1) 1 1 
i) Qk—-1 F245 114) ’ ’ _ 2), 
ke? (5 — pi) + olla — 5) (7) 
2k+1 
For k + oo, making the approximation that ( . wf EPR re a Vk 


F,(2,1—k;k+3;-1 
and 2/1 a dn ‘, then 


1 k 1 
Pen) os aa »fid — pi). (8) 


. QPyC 
That means for given Paces 


2 

ae ae 

mF ( oe ©) 
2 


yielding a critical exponent of 2. 

We will now compare the performance of QPyC with surface codes [24] for the 
correction of erasure errors. Surface codes consists of a D x D square lattice with 2D? 
qubits in total, where each qubit located on the edges of the square lattice. The logical 
operators X;, (Z,) of the surface codes are given by the product of X (Z) operators 
along a non-trivial homological cycle connecting the boundaries [24]. Therefore, the 
success probability of bond percolation in the square lattice and its dual lattice [22] is 
essentially the same as the probability of decoding the surface code. It is well known 
that the bond percolation threshold is 50% for a square lattice, meaning that as the 
distance of the lattice D —> 00, P®°..., > 1 for p, < 50% and PS... + 0 for p > 50%. 

In Fig. [2(b)] we study the success probabilities of surface codes with three different 
distances D = (3,7,11) and find that at p, = 50%, P&C... © 0.30, while in Fig. |2(a) 


success 


we see that P2EYC. = 0.5 for all code sizes. Therefore, when p; is around 50%, QPyC 
of any code size always outperforms surface codes. Further, it can be seen in Fig. 
that threshold for surface codes with D = 3 is about 37%, and the threshold approaches 
50% for larger code sizes, while the threshold of QPyC is 50% for all code sizes. In 
practice, we might operate the codes at a loss probability much lower than 50%. Hence, 
we assume p; © 20% and compare the surface code and QPyC with 22° ~ (2k +1)24+) 
so that they have a similar Hilbert space dimension for physical encoding. As illustrated 
in Table 1, for smaller code sizes surface codes performs better than QPyC and with 
D>9andk > 15, QPyC outperforms surface codes with a smaller failure probability. 


4. Application in quantum repeaters 


One-way (or third generation) quantum repeaters rely on quantum 
error correction to relay data from one repeater station to the next. At each station, 
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k; | Surface code QPyC 

6 0.0068 0.007 

oar x Ii 0.0016 

15| 12x10 8.8 x 1075 
21) 12410" | 523% 10° 


monoy 
o 


Table 1. A comparison between the failure probabilities of error correction (1 — 
Psaccess) for surface codes and QPyC with a similar size of the Hilbert space for 20% 
erasure rates. All values were calculated for surface codes with a total of 10° runs. 


Figure 2. (a) Success probability of QPyC for different erasure rates. (b) Success 
probability of decoding surface codes with three different distances (3,7,11). 


error correction operations are performed before the message is transmitted to the 
next station. Recently, it has been shown that QPC (a generalization 
of the |[4, 1, 2]] code) can correct erasure errors efficiently by teleportation based error 
correction (TEC). Since QPyC corrects a larger fraction of erasure errors, it is reasonable 
to consider QPyC instead of QPC for one-way quantum repeaters, and expect a 
significant improvement for resource consumption and key generation rates. 


4.1. Implementation of error correction for erasure errors 


For qubit encoding schemes, it has been shown that TEC [83} [34] is an effective approach 
for the correction of erasure and operation errors. Since teleportation requires X and Z 
measurements on qubits, error correction for erasure and operation errors is possible if 
we can reliably measure the logical operators X;, and Z_ of the error correcting code [26]. 
For example, the |[4,2,2]] code has stabilizers XX XX and ZZZZ and the following 
logical operators: X{) = IXIX,XIXI, XP = 1Z1Z,Z1ZI, Z” = ZZII,IIZZ and 
Zo) = IIXX,XXII. If any qubit undergoes an erasure error, all the X; operators 
and Zz, operators can be measured in the TEC circuit and the encoded qubits can 
be retrieved. For example, if the first qubit is lost, one can calculate the operators 
X = 1X1X, X° = 1212, Z = 11ZZ and Z = IIXX and retrive the encoded 
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Figure 3. TEC circuit for multilevel systems, where erasure and operation (X and Z) 
errors in |q) 7 are corrected. The states |0); and |+); = pa |7) should be prepared 
fault-tolerantly and free from erasure errors. ”S” refers to an encoded SUM gate that 
acts on logical qudits |¢); and |j)z, as |¢)r\j)z — |¢)x|(¢ + 7)mod d);. The encoded 
SUM gate has pairwise implementation for CSS codes. 


qubits. We generalize TEC to multilevel systems in Fig. by using generalized 
Pauli matrices that act on d-level system as X'|j) = |j + 1) and Z'|j) = w)|j), 
0<ij<d-1 and SUM gate that acts on a control qudit |/) with a target 
qudit |j) to produce the transformation |7)|j) — |i)|(i +7) mod d) [85]. Consider 
the [[3,1,2]]; code with stabilizers XXX, ZZZ and X_ operator that corresponds to 
one of the operators |X X?, X XI, X?IX and Z, operator that corresponds to one of 
the operators [Z?Z, ZIZ?, Z*7ZI. It can be readily seen that in the presence of an 
erasure error on any qutrit one can still measure the corresponding logical operator and 
reconstruct the encoded state. 


4.2. Error model for operation errors 


Suppose that each physical qudit is encoded into a [[2k + 1,1,k + l]]a QPyC and is 
transmitted through repeater stations, where the TEC is implemented. We assume 
independent errors acting on physical qudits and consider an extensive error model 
with the following types of errors acting on the encoded qudit as follows: 


eee 
(i) The photon arrives at each repeater station with probability 1—p, =e “att , where 


Lo is the repeater spacing and Lay_ = 20km is the attenuation length of the fiber. 

(ii) The photon undergoes depolarization with a probability eg. For simplicity, we 
assume that this is also the probability of error on each physical qudit of the encoded 
Bell pair prepared at the repeater station. 

(iii) The photon experiences an additional dephasing in the matter qubit-photon 
coupling with a probability €,. 

(iv) Each SUM gate acting between the encoded Bell pair and the incoming qubit (see 
Fig. |3) fails at a probability €,. 

By further assuming the same probabilities for all types of depolarization errors and 

gate errors, the transmission channel with the incoming single qudit state A takes the 


Overcoming erasure errors with multilevel systems 9 


(10) 


form, 
(l- peo 
E.(pa) = (1 — pi)(1 — €a ~ &)pa + pilvac) (vac| + —— : (X'Z!)pa(X'2?)' 
1,7=0 
an) : d-1 
l rs Ze \pal (z*y! 
k=0 


Similarly, the imperfect gate between the incoming qudit A and qudit B at repeater 
station can be modeled as 
Eq(pas) = (1 ~~ en UP aB Ue, 
d-1 
€ 
+ S07 (X°2X2*) pap(XAZIXAZ4), (11) 
a,b,c,e=0 
For each qudit the locally prepared logical states |0); and |+), = ae lj), we 
assume it undergoes depolarization, which can be modelled as 
d-1 
E,(pp) = (1— «apn + $ Zp = > (X12) pp(X'Z7)'. (12) 
1,j=0 
In contrast to the error models in [10], we track the errors in the preparation of 
the logical Bell state needed for TEC. From these error channels, we can calculate the 
probability of having an error in os one of the measurements up to the first order as 


SE 
€x = race — d*) + @ ie ee d) 
3€ € 
€z = aE (d* — d®) + “ eg d) + aid —1). (13) 


Once the encoded state [[2k + 1,1,k + 1]]q is transmitted from one repeater station to 
its neighbor, there exists three possibilities at the receiving repeater station 


(i) More than & photons are lost in transit and the outcome of the measurement cannot 
be found leading to a heralded failure with probability Pai = 1 — PQPyC | 


(ii) At least & photons are received, but the encoded state is not decoded correctly 
due to the presence of many operation errors and the encoded state is not decoded 
correctly with probability Pincorrect: 


(iii) At least k photons are received to make an encoded X/Z measurement and the 
encoded state is decoded correctly with probability Pooprect: 


The probability that a heralded failure does not happen at any one of the r repeater 
stations is given by, IPoe onl . Let us suppose that n, photons are lost before the 


Success 


destination. Among the rest of the (2k + 1 —n,) photons that reach the destination, n2 
photons suffer operation errors. As such, the code can correct up to nj +2n2 < k errors. 
The probability of successfully measuring the encoded X/Z measurement outcomes is 


given by, 
2k+1\ (/2k+1—-—n 
Preorrect( (X/Z) = 3 csi fos ) ( ~ ) x 


ny= =0 n2= =0 
2k+1-— a al 2k+1—n1—n2 (14) 


(pi) €X77,(1 — pr) — €x/z) 
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The probability of incorrect decoding of the qudit is given by, 


2k+1—n, 
2k+1 2k+1- N41 
Pincorrect( (X/Z) = 3 > ( i ) ( - ) x 


m=Ong=[5- F451 


\2ke1=m (1 = exjz) tm (15) 


(pi) x77, — pt 
It is easy to verify that Praicx/z) + Poorrect(x/z) + Pincorrect(x/z) = 1. By making a 
pessimistic assumption that an effective logical error in any one of the repeater stations 
is an overall logical error at the receiver’s end, the logical error rate of the encoded 
quantum information can be defined conditioned on the success of receving enough 
photons as, 
[ Poorrect(X/Z) | 
[Paes 


Success 


Qx/z=1- (16) 


For the two basis protocol for quantum key distribution (where information is encoded 
in only two logical bases X! and Z'), the asymptotic normalized secure key generation 


rate is [36]. 


(Poe 
i= oo (logod — 2h(Q)), (17) 
where to is the time taken for local operations with 
@- 2 Ge 5 22) 
h(Q) = —Qlogy 7 ~ (1 — Q)logy(1— Q). (is) 


In Fig. 4(a), we study the dependence of the secure key generation rate Rto, 
with respect to the total distance of communication for small encoded blocks of QPyC, 
assuming 1 km repeater spacing and no operation errors. For to = 1 ps, 3 qudits are 
sufficient to reach 700 km with key generation rates R ~ 10kHz. We can increase 
the total range of communication to 10,000km by using 7 qudits and obtain a key 
generation rate R = 1000kHz. As the operation errors involved with the error correction 
€g = €d = & = € increases, Pincorrect(x/z) Increases exponentially after a certain distance 
Diop and we can expect a quick decay of R.tg over a certain Li. The maximum Q 
that the QPyC can tolerate depends on the dimension d of the code. We can check 
that Qmax 0.15, 0.21, 0.237 for the [[3, 1, 2]]3 code, [[5,1,3]]; code and ||7, 1, 4]]7 codes 
respectively for the codes to yield non zero key generation rates. Taking only the leading 
terms into account in Eq. (14), we can approximate Prncorrect(X/2) & 3ex/z(1 — pi) 
for the [[3,1,2]]3 code, Pincorrect(X/2) © dex zp} \(a —p,)* for the [[5,1,3]]; code and 
Prncorrect(X/Z) © 5ex/z(pi)? (4) (1 —p,)° for the {[7,1, 4]]7 code respectively. Further, since 
Pra) for 1km spacing is negligible, we should have Pineorrect X 2% = Qmax, Yielding 


Lo 
L a 2QmaxLo 
tot,max ™ 
Fcoreeci( = F rpereet a) 


(19) 


For example, for « = 10°-*, Lo = 1km, we get Litotmax = 120km, 440km, 1900km for 
[[3, 1, 2]}3, [[5, 1,3]]; and [[7,1,4]]7 codes respectively as confirmed by the ” x” lines in 
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Fig. (4). In Fig. 4(b-d), we study the variation of R.to for varying operation errors 
Eg = €4 = & = 10°°, 10-°, 10 for [[7, 1, 4]]z, [[5, 1, 3]Js and [[3, 1, 2]]; codes respectively. 
We fit the points obtained from a rigorous theoretical calculation with the approximate 
Pincorrect Mentioned above (shown as gray curve in Fig. 4) and observe a good match 
between them. 
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Figure 4. (a) The key generation rates Rto that can be achieved with small blocks 
of QPyC in the absence of operation errors with 1km repeater spacing between the 
repeater stations with [[7,1,4]]7 (red), [[5,1,3]]5 (blue) and [[8,1,2]]3 (green) codes. 
(b-d) R.to in the presence of operation errors €, = €¢ = €» = € = 10°, 10~°, 10~* for 
([7,1,4]]7 (b), [[8,1, 2]]s (c) and [[5,1,3]]5 (d) codes. The gray line corresponds to the 
approximation taking only the leading terms into account. 


4.3. Comparison with QPC 


We can compare the performance of different quantum codes for one-way quantum 
repeaters by considering both the qubit and temporal resources consumed by the code, 
respectively. In order to compare QPyC codes, which use multilevel systems, with QPC 
of qubits, we consider the conversion of a d-level qudit into [log,d] qubits and compare 
the performance of the codes by a cost coefficient 


i 2(2k + 1)/log,d| 
ae LoR 


qubits/km/sbit/s, (20) 


The Cj, is obtained from the product of qubit resources and temporal resources. The 
number of qubits used for TEC at every repeater station is given by 2(2k + 1)|log,d| 
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and the number of qubits for all stations is 2(2k + 1)[logod](Ltot/L0) . The temporal 
resource used by the scheme is simply the inverse of key generation rate, i.e. (1/R). 
Since, the product of qubit and temporal resources varies at least linearly with Li.., we 
further divide the product by Lj, to obtain the cost coefficient, which stands for the 
number of qubits required per km for the generation of one secure bit in one second. 
We find that the present scheme with QPyC can achieve a very small cost coefficient, 
which is about 5 times less than for QPC with TEC for £4, = 10,000 km in the absence 
of operation errors. Note that in addition to the local resource overhead, we can also 
save the number of modes considerably by using QPyC instead of QPC (See Appendix 
A). In the presence of operation errors, the comparison between qubit and qudit based 
schemes for QRs depends largely on the error model because of the complexity involved 
in the implementation of multimode operations. Although it is challenging to compare 
the resource requirements between these codes, we have attempted to ensure a fair 
comparison by assuming that the operation errors increase with dimension d of the 
qudit. We assume that €, = €, x d*, eq = €j x d’, and €, = & x d for the comparison. 
In Fig. |5| we investigate the variation of the ratio of cost coefficients Copg/Copyc with 
respect to total distance of communication L;,;. We consider only one of the errors €,, 
€q and €, respectively in each one of the plots and show the variation of the ratio of the 
cost coefficients with respect to L;.; and the corresponding d picked by the optimization 
of the cost coefficient shown in the colored contour. As the operation error increases, 
the ratio decreases and QPC becomes more favorable than QPyC shown in the area 
above the break-even contour line (shown as thickened line) as the number of levels d 
increases. Since €, scales as d* and d > 2k +1, the optimization is forced to choose a 
smaller code for large total distances in Fig. 5(a) compared to Fig. 5(b) and Fig. 5(c). 


5. Physical implementation 


For the implementation of quantum repeaters with QPyC encoding schemes, we can 
consider time-bin photonic qudits. The time-bin qudit has a single photon excitation in 
a superposition of d time-bins, which can be efficiently generated at telecom wavelengths 
and coupled into an optical fiber [37]. The TEC circuit, which is essential of our scheme, 
requires single qudit X- and Z-operation, two-qudit SUM gate, X- and Z-measurement 
of qudits. The SUM gate can be further decomposed using Controlled Z (CZ) gate and 
Fourier gates [35]. 

For the time-bin qudits, the Z-operation and Z-measurement can be implemented 
by a selectively phase shifting or detection at different time slots. X-operation and 
X-measurement can be achieved by making a strong dispersive medium that interferes 
different time-bins as detailed below. For the two-qudit gate, the direct photon-photon 
interaction is negligible due to weak single photon nonlinearity. Therefore, we propose 
to use a nanophotonic-atom interference to store and manipulate the photonic qudits 
and to mediate photon-photon qudits gate. 

The feasible experimental system adapted from Ref. [88] is illustrated in Fig. 6(a). 
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Figure 5. Contour plots (a), (b) and (c) show the improvement factor of QPyC 
with respect to QPC Copc/Capyc when the imperfections are dominated by the two- 
qudit gate error €, = €,/ d*, qudit depolarization error €y = €a/d? and qudit dephasing 
error €) = €,/d respectively. The optimal d chosen by the optimization (of C’) is 
shown with different colors in the same plots. The contour area below the break even 
contour line (shown with thicker line) indicates the area where QPyC performs better 
than QPC. Contour plots (d-f) show the variation of optimized cost coefficients Cj, for 
different dominant errors. Here, it is assumed that it takes the same time to create 
small encoded block of qubits(qudits). 


A single atom is trapped in the vicinity of a fiber integrated photonic crystal nanocavity, 
single photons can be efficiently coupled to the nanocavity using the tapered fiber, which 
strongly interacts with the atom due to the very small mode volume of the nanocavity. 

Based on the nanophotonic atom interface, we can implement the Fourier (F’) 
and X-gates of single time-bin qudit by downloading the photonic time-bin qudit to 
the atom, and using photonic-atomic qudit control phase (CZ) gate. Based on these 
elementary gates, the X-measurement can be realized by the combination of F-gate 
and Z-measurement. The photon-photon qudit SUM gate can be realized by a sequence 
of Fourier gate on photon qudits and two CZ gates between photon and atom. In 
the following, we provide the details about experimental implementation of the above 
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element gates: (1) F and X gate of single qudit, (2) the photon-atom qudit CZ gate, 
(3) the atom mediated photon-photon gates between time-bin qudits. Based on these 
three building blocks, QPyC can be corrected for erasure and operation errors using 
TEC and retransmitted to the neighboring station. 

To illustrate the procedure, we consider the simplest case of qutrits. We choose 
the Dy lines of natural ®’Rb atoms [39], with the transitions between hyperfine levels 
|So,1,2) = |5’Sij2, F = 1, mp = —1,0,1) and |5?P32, F = 2, mp =0,1,2). In addition, 
we use the |A) = |57S,/2, F = 2, mp = —2) as an ancillary state for processing the 
photonic qutrit. As illustrated in Fig. 6(b), we apply a tunable external magnetic field 
to introduce a relative frequency shift depending on hyperfine levels, which will be useful 
for the atom-photon coupling gate. 


5.1. Single qudit F and X gates 


Since it is difficult to make an X-measurement on a time-bin photonic qutrit, we propose 
to transfer the quantum state from the time-bin photonic qutrit to either a multilevel 
atom or multiple two-level atoms. First of all, with microwave driving of transitions 
between F=1 and F=2 levels and also optical pumping, the single °’Rb atom can be 
initialized to the ancillary state. Then, for single photon input to the cavity, it will 
be largely detuned from the transitions of the atoms. By applying an external laser 
driving with frequency wr — w, = Apr, where Ap is the energy difference between |A) 
and |.S;), through a stimulated Raman adiabatic passage [40], the population will be 
transferred to |S,) if there is a single photon in the incoming time-bin, as shown in 
Fig. 6(b). By a sequence of microwave pulses, the population will be mapped to |.5;) 
according to the order of input time-bin. This process effectively maps the state of 
the time-bin qutrit to the atom’s ground state levels. After that, we can carry out an 
X-measurement by performing a Fourier transformation followed by a Z measurement 
on the atom. Alternatively, it is also possible to store the time-bin qutrit into multiple 
atoms and control the process by synchronized pulse sequences on individual atoms. 

This procedure can be generalized to time-bin qudits (d > 3) in a straightforward 
fashion. When the time-bin qudits are stored in atoms, arbitrary single qudit operations 
can be realized by applying either optical or microwave pulses. In the most general case, 
arbitrary unitary on single d—level system can be constructed by a series of SU(2) gates 
between any two of d levels (called Givens rotations) [41]. Alternatively, the arbitrary 
unitary can be realized with optimal controlled pulse sequences. For the qutrit case, 
there are more efficient schemes that only requires three steps, by either three Givens 
rotation or three Houlsehould’s reflections [42]. Recently, arbitrary unitary on d = 16 
hyperfine energy levels of 13°C's atom has been demonstrated [43]. Two-qudit gates 
between atoms can be mediated by cavity photons. Generalization of such techniques 
will enable the generation of the encoded Bell states needed for the TEC at repeater 
stations. 
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Figure 6. a) The schematic of a single atom trapped close to a photonic crystal 
nanocavity and the time-bin photonic qutrit input from the fiber interacting with the 
atom. b) The energy diagram of transfer of the quantum state of a time-bin qutrit to 
an atom. The atom state is initialized to the |A), the input photon and the strong 
driving (Red arrow) are detuned from the excited state by Ar. Laser pulse followed 
by a joint microwave pulse is applied to map the state of the time-bin qutrit to the 
atom. c) The energy diagram of the CZ gate between photonic qutrit and atom. 


5.2. Controlled-Z gate between photonic and atomic qudits 


To realize CZ gate between the time-bin photonic qutrit and the atom, we prepare the 
atom in the superposition of |S;). Due to the strong cooperative interaction between 
the atom and nanocavity photon, the incoming single photon pulse input to the system 


have the reflectivity as 
=e 
pt ey (21) 
n+1— 2id/y 


where 7 is the atom-photon interaction cooperativity, 6 is the frequency detuning 


between the pulse and atomic transition frequency and y is the excited state decay 
rate. For a large n(~ 100), r ~ e*? and @ & 2arctan =e Due to the magnetic field, 
the photon detuning 4; for different |.S;) gradually reduces with j, corresponding to a 
atomic state dependent phase shift as shown in Fig. 6(c). 

The CZ gate can be represented as C = 5°, |m)(m| ® Uam, where m = 0,1,2 
denoting the states of photonic qutrits and U,,, is photon state dependent phase gate. 


The gate can be realized by the following three steps: 


(i) For the first incoming pulse, all energy levels are detuned from the cavity mode 
by controlling external bias magnetic field such that d9 = 6; = 62 then we have an 
operation on the atomic energy levels as Ua = diag{1, 1, 1}. 
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(ii) For the second pulse, magnetic field is tuned to have 69 = V3yn/2, 6, = 0 and 62 = 
—/3yn/2, so that the phase shift is 27/3, 0 and —27/3, Ug. = diag{e's ,1,e7*5 }. 

(iii) For the third pulse, the magnetic field is reversed, so that the phase shift is —27/3, 
0 and 27/3, then Ua2 = diag{e-*3 , 1, e's }. 


The generalization of the CZ gate to higher dimensional qudits (d > 3) is less 
straightforward compared to time-bin storage. One possible realization is using the 
cavity coupled multiple A-type atoms (each atom is an effective two-level system). 
Suppose there are d atoms and each has ground states |g), |s) and excited state |e), the 
transition |g) < |e) is near resonance with cavity while |s) < |e) is far off-resonance. 
In addition, the transition frequencies of atoms can be controlled by external electric or 
magnetic field individually. First, the atoms are initialized to |ss...s), and transition 
frequencies are on-resonance with cavity. By single photon Raman transition of the 
atom ensemble, the states of atoms are prepared to the one-excitation Dicke state 
allgs- ..8)+|sg...s)+...|ss...g)). Then, the CZ gate of d-level time-bin qudit can be 
realized by shifting the transition frequencies of individual atoms: for i-th (i = 1,...,d) 
pulse input to the system, the j-th (j = 1,...,d) atom is tuned to be near-resonance 
on the cavity so that the reflected pulse gains a phase of e~7a(—DU-D, 


5.3. Controlled-Z gate between time-bin qudits 


A CZ gate can be realized between two physical time-bin photonic qudits by a 
generalization of the Duan-Kimble scheme [7]. The CZ gate between two time 
bin photonic qutrits (f and s) can be realized by preparing the atomic state in a 
equal superposition of all energy levels and applying three photon-atom CZ gates 
(Pes = dom lM) p5(M| ® Uam with m = 0,1,2 denoting the states of photonic qutrits) 
and two Fourier gates (F’) on atoms as 


Uog = Cr °F *C, °F G;, (22) 
where, 
a 1 
fs A 1 ee ete (23) 
1 es es 


The atomic state initialization and Fourier gates of atoms can be realized by microwave 
or optical Raman pulsed controlled transitions (|So) < |.S1) and |S,) < |S2)) and phase 
gate Z. 

For qudits with d > 3, the Fourier gate of atoms can be realized by virtual 
cavity photon mediated atom-atom interaction, where strong pumping is applied on 
the |s) + |e) transition, with the detuning equals to that between cavity photon and 
|g) < |e). Based on this procedure and the CZ gate between photonic and atomic qudits 
(d > 3) described in the previous section, the CZ gate between two time-bin qudits can 
be realized using a similar technique as Eq. (22). 
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6. Conclusion 


We have investigated efficient codes using multilevel systems that can correct up to 50% 
erasure error rates, which is the bound set by the no-cloning theorem [I]. The success 
probability of quantum polynomial codes close to 50% erasure rates is higher than 
the success probability of surface codes of all code sizes. We employed teleportation 
based error correction that can correct erasure errors up to the threshold efficiently 
and discussed its application for the construction of highly efficient one-way quantum 
repeater networks. In comparison with quantum parity codes, we have obtained an 
improvement in resources by about 5 times for communication across 10,000km in the 
absence of operation errors and numerically identified the parameter regime, where the 
quantum polynomial codes perform better than quantum parity codes in the presence of 
operation errors. We have discussed the physical implementation of quantum polynomial 
codes and identified the key technological requirements. It will also be interesting to 
consider the experimental implementation of multilevel systems using oscillators [44], 
Rydberg atoms and ensembles of multilevel systems [46]. Moreover, we can extend 
the coding schemes to continuous variables for the correction of erasure errors [47] [48]. 
Besides quantum communication, quantum error correcting codes for multilevel systems 
that can correct a large fraction of erasure errors might be useful for improving precision 


metrology [19] BS]. 
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APPENDIX A 


In the manuscript, we have assumed that we are only limited by the qubit resources 
rather than the optical modes of the fiber channel. If qubit resource is no longer the 
limiting factor for one-way quantum repeaters, we should then compare the number of 
modes needed for different repeater schemes. For quantum state transfer using time- 
bin qudits, each qudit needs d temporal modes for transmission. Therefore, the cost 
coefficient based on the number of modes for the transmission of photons is given by, 


' __ (2k+1)d 
Crn(Ltot) = fees eS 


In Fig. |7} we compare the cost coefficients Ci, of QPyC and QPC. For QPC, since each 
qubit needs two temporal modes for transmission the numerical values of C, and C!, are 


modes/km/sbit /s. (24) 
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the same. For QPyC, the numerical value of C’, is increased by a small factor compared 
to Go. We find that for L4., = 10,000 km, QPyC achieves a cost that is about 3 times 
less than QPC. 
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Figure 7. The cost coefficient C/, based on the number of modes for QPC and QPyC. 
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